mat  wdhcimyl  ..-..•> 

AFCRC-TN-60- 1  1  1 

25  Waverly  PI**,  New  York  J,  N.  % 


£\     y-r    A  NEW    YORK    UNIVERSITY 

Su        |sJ  ^  Institute  of  Mathematical  Sciences 

^      V-y        <  Division  of  Electromagnetic  Research 

^CCCXX1" 


RESEARCH    REPORT    No.    EM-150 


Asymptotic  Solution  of  Systems  of 
Linear  Ordinary  Differential 
Equations  with  Discontinuous  Coefficients 


CLIVE  R.  CHESTER  and  JOSEPH  B.  KELLER 


Contract    No.    AF    19(604)5238 
DECEMBER,    1959 


AFCRC-TN- 60-111 


NEW  YORK  UNIVERSITY 

Institute  of  Mathematical  Sciences 

Division  of  Electromagnetic  Research 


Research  Report  No.  EM-150 


ASYMPTOTIC  SOLUTION  OF  SYSTEMS  OF  LINEAR  ORDINARY 

DIFFERENTIAL  EQUATIONS  WITH  DISCONTINUOUS 

COEFFICIENTS 

Clive  R.  Chester  and  Joseph  B.  Keller 


Clive  R.  Chester 


Joseph  B  Keller 


Morris  Kline 
Director 


MnuW. 


Dr.  Werner  Ger 
Contract  Monitor 


¥* 


The  research  reported  in  this  document  has  been  sponsored 
by  the  Electronics  Research  Directorate  of  the  Air  Force 
Cambridge  Research  Center,  Air  Research  and  Development 
Command,  under  Contract  No.  AF  19(6o4)52j58. 


Requests  for  additional  copies  by  Agencies  of  the  Department  of 

Defense,  their  contractors,  and  other  Government  agencies  should 

be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
DOCUMENTS  SERVICE  CENTER 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA 

services  or  have  their  ' need-to-know '  certified  by  the  cognizant 

military  agency  of  their  project  or  contract. 

All  other  persons  and  organizations  should  apply  to  the: 

U.S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHINGTON  25,  D.C. 


Abstract 

A  system  of  linear  first  order  ordinary  differential  equations  is 
considered  in  which  the  coefficient  matrix  depends  upon  a  parameter  k. 
The  asymptotic  form  of  the  solutions  of  this  system  are  determined  for 
large  values  of  k  when  the  coefficient  matrix  or  some  of  its  derivatives 
have  finite  jump  discontinuities.   The  asymptotic  form  obtained  contains 
terms  beyond  those  previously  known.   These  new  terms  have  discontinuous 
coefficients,  whereas  the  previously  known  terms  are  continuous. 

The  results  are  applied  to  the  one  dimensional  propagation  of  time- 
harmonic  waves.  It  is  shown  that  in  this  case  the  reflection  coefficient 
is  proportional  to  k"   if  the  J-th  derivative  of  the  index  of  refraction 
is  discontinuous  and  the  lower  order  derivates  are  continuous.   Here 
k  =  2ir/x  is  the  propagation  constant  and  X  the  wavelength  in  a  medium  of 
unit  refractive  index.   This  result  is  to  be  contrasted  with  the  fact 
that  the  reflection  coefficient  vanishes  faster  than  any  negative  power 
of  k  if  all  derivatives  of  the  index  of  refraction  are  continuous. 
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1.   Introduction 

Let  us  consider  the  vector  solution  u(x,k)  of  the  system  of  N  first 
order  linear  ordinary  differential  equations 

(1)         g  =  kA(x,k)u 

In  this  system  the  coefficient  matrix  kA  depends  upon  a  real  parameter  k. 
Asymptotic  expansions  of  u  in  powers  of  k   ,  which  are  valid  as  k  -*  oo  , 
are  well  known.   They  hold  provided  that  A  has  continuous  derivatives  of  all 
orders  and  satisfies  certain  other  conditions.   But  if  some  derivative  of  A 
is  discontinuous,  then  only  a  finite  number  of  terms  in  the  expansion  of  u  are 
known.  We  consider  the  problem  of  finding  additional  terms  when  A  and  its 
derivatives  have  finite  jump  discontinuities.  We  shall  show  how  to  solve  this 
problem,  and  we  shall  find  explicitly  the  next  term  beyond  those  previously 
known.   This  term  is  discontinuous,  whereas  the  preceding  terms  are  all  con- 
tinuous.  This  fact  is  significant  in  various  applications. 

One  important  application  of  our  result  is  to  one  dimensional  propagation 
of  time  harmonic  waves  in  an  inhomogeneous  medium.   Such  a  wave  motion  w  satis- 
fies the  reduced  wave  equation 


(2)       d2w   ,2  2,  , 

v  '       — p  +  k  n  (x)w  =  0. 

dx 


Here  n(x)  denotes  the  refractive  index  of  the  medium  and  k  denotes  the  propa- 
gation constant  in  a  medium  with  n  =  1.   It  follows  from  known  results  that 
the  reflection  coefficient  R(k)  determined  from  (2)  is  asymptotically  zero  to 
all  orders  in  k   ,  provided  that  n  has  continuous  derivatives  of  all  orders. 
However,  if  the  J-th  derivative  of  n  has  a  finite  jump  discontinuity,  and  the 
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lower  order  derivatives  are  continuous,  our  result  shows  that  R(k)  is  of  the 
order  k   .   This  essential  difference  between  the  infinitely  and  finitely  dif- 
ferentiable  cases  has  been  observed  by  S.A.  Schelkunoff  *-  ^  and  J.  Feinstein'-  -*  • 
Equation  (2)  is  treated  in  section  h. 


2.    Preliminaries 

We  wish  to  determine  the  asymptotic  behavior,  for  large  values  of  k,  of 
the  solutions  of  (l)  in  an  interval  I  of  the  real  x  axis.   We  assume  that  with- 
in I  the  matrix  A  has  the  asymptotic  expansion 


(3)       A(x,k)  -^  Aj(x)k" 


J=o 


We  also  assume  that  A  (x)  has  N  eigenvalues  X  (x) , . . . X^(x)  which  are  distinct 
in  I.   Then  there  exists  a  non  singular  matrix  R(x)  which  diagonalizes  A  (x) 
in  I.   We  define  the  diagonal  matrix  A(x)  by  the  equation 

(!)■)       A(x)  =R"1AoRsdiag[l1(x),...,yx)]. 

The  matrix  R  can  be  so  chosen  that,  like  A,  it  is  independent  of  k  and  possesses 
as  many  continuous  derivatives  as  does  A  (x) . 

With  the  aid  of  R  we  define  the  vector  v(x,k)  by  the  equation 

(5)       u(x,k)  =  R(x)v(x,k). 
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Then  from  (l),  it  follows  that  v.  satisfies  the  equation 


(6)       ^=kB(x,k)v. 


Here  the  matrix  B  is   defined  by 


(7) 


B(x,k)    =  R_1AR-k  ■LR_1dR/dx. 


It  has  the  asymptotic  expansion 


(8) 


}(x,k)  ~  )   B.(x)k  J. 


J=o 


The  B .  are  given  by 

(9)       B  =  A,  B  =  R_1A^R-R"1dR/dx,  B  =  R_1A.R,  j  =  2,3, ..  . 

We  shall  determine  the  asymptotic  behavior  of  the  solutions  v  of  (6).   Then 
the  behavior  of  the  solutions  u  of  (l)  will  be  found  by  means  of  (5).   The 
advantage  of  dealing  with  (6)  rather  than  (l)  is  that  B  is  diagonal  while 
A  is  not  necessarily  diagonal. 

Let  us  tentatively  assume  that  (6)  has  N  linearly  independent  solutions 
v.(x,k)  with  asymptotic  expansions  of  the  form 


(10) 


v±(x,k)  ~  exp 


k   /   Xi(x)dx  I    \   v°(x)k"S]  i=l,...,N. 


x 

*-    o 


s=o 


Let  us  also  suppose  that  the  expansion  of  dv/dx  is  obtained  by  termwise  dif- 
ferentiation of  (10).   Then  we  can  insert  (10)  and  (8)  into  (6)  and  equate  to 
zero  the  coefficient  of  each  power  of  k   .   In  this  way  we  pbtain  the  following 
set  of  equations: 
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s-x 

(11)      M.vS  =  dvS_1/dx  -  >    B   .vJ 
11     i        /_,    s-j  i 


B-l 

i 

s  =  0,1, . . . 


Quantities  with  negative  subscripts  or  superscripts  in  (ll)  and  elsewhere  are 
defined  to  be  zero,,  and  M.  is  defined  by 

(12)  M.  =  Bq-  X.I  =  diag^-  X.,...,^-  X.]. 

M.  is  diagonal  with  a  zero  in  the  (i,i)  position  and  it  is  therefore  a  singular 
matrix. 

If  we  set  s  =  0  in  (ll)  we  see  that  the  right  side  vanishes.   Therefore 
the  only  non  zero  component  of  v.  is  the  i   one.   Then  the  solvability  condi- 
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tion  for  (ll)  with  s  j=   0  is  simply  that  the  i   component  of  the  right  side 
vanish.   Upon  replacing  s  by  s  +  1,  this  condition  can  be  written  in  the  form 

s  1 

(13)  dv^./dx  -  (B1v^)i  =  ^(Bs+1_(jvj).      s  =  0,1,-..  • 

The  second  subscript  indicates  the  component  of  the  vector.  Equation  (13)  is 

a  first  order  ordinary  differential  equation  for  v. . .   From  (ll)  all  other  com- 

s  1s 

ponents  of  v.  can  be  determined  in  terms  of  the  vV  with  j  <  s.   Then  v^  can  be 

found  from  (13) .   Thus  by  starting  with  s  =  0,  we  can  determine  successively 

s 
all  the  v. . 

1 

By  keeping  track  of  the  differentiation  which  is  performed  in  (ll)  at  each 
step  of  this  process,  we  can  determine  the  continuity  properties  of  the  deriva- 
tives of  the  v?  in  terms  of  those  of  the  B.  and  therefore  of  the  A  .   In  this 

J  J  r  -1 

way  we  arrive  at  the  following  lemma,  which  is  due  to  P.D.  Tamarkin  L  J . 


Lemma  If  for  some  J  *  0  and.  all  j  £  J  ,  each  coef- 
o  o 

ficient  A  has  J  -  j  continuous  derivatives,  then 
J      ° 

so  does  v.  for  i  »  1, ...,N. 

This  lemma  can  be  verified  readily  by  induction,  with  the  aid  of  (ll)  and  (lj) 
Tamarkin  has  also  proved  the  following  theorem. 

Theorem  I   Suppose  that  in  an  interval  I  of  the  x  axis 
A(x,k)  has  the  asymptotic  form 
J 


(Ik)  A(x,k)  =  V  A.(x)K~J+   o(k  °) 


Let  the  A.  satisfy  the  hypothesis  of  the  lemma  and  let 

A  (x)  have  distinct  eigenvalues  \.(x),    i  =  1,...,N  in 

I.   Let  uf(x)  =  Rv.-  where  v.  are  the  vectors  constructed 
j.x  i        1 

as  described  above  and  R  is  a  matrix  which  diagonalizes 

A  and  has  as  many  continuous  derivatives  as  does  A  . 
o  ° 

Then  there  exist  N  linearly  independent  solutions  of 

(l),  u.(x,k),  i  =  1, ...,N,  such  that  in  I 

J 
x  o  -J  \ 

(15)      u.(x;k)  =  exp  k  /  Xi(x)dx  f  )      u^(x)k"S+  o(k  °)J 


j=o 


s 
The  coefficient  u  (x)  has  J  -  s  continuous  derivatives. 

3.   The  discontinuous  case 

Let  us  now  suppose  that  J  is  the  largest  value  of  J   for  which  the  hypo- 
thesis of  the  lemma  is  satisfied.   Then  for  at  least  one  value  of  j  in  the 
range  0  &    j  §  J  +  1,  the  coefficient  A.  has  a  discontinuous  (J  +  1  -  j)th 

J 


derivative.   Let  us  suppose  that  this  discontinuity  occurs  at  x  .   Then  ve  may 
choose  I  so  small  that  for  each  j  in  the  range  0  §  j  §  J  +  1,  the  (j  +  1  -  j)-th 
derivative  of  A.  is  continuous  in  I,  except  perhaps  at  x  .   We  further  assume 
that  the  only  discontinuities  of  any  of  these  derivatives  are  finite  jump  discon- 
tinuities at  x  .   From  the  definition  (9)  of  the  B.  if  follows  that  the  proper- 
ties just  described  pertain  to  them  also. 

Let  us  now  consider  the  solutions  v(x,k)  of  (6)  in  which  the  coefficient 
matrix  is  B(x,k) .   With  the  asumption  just  made  on  the  B.,  the  hypotheses  of 
the  lemma  and  of  Theorem  I  are  satisfied  by  B  with  J  =  J  +  1,  in  each  of  the 

two  subintervals  I   and  I   into  which  x  divides  I.   Therefore  there  are  N 

o 

linearly  independent  solutions  v.  of  (6)  which  have  the  asymptotic  form  (15)  in 

I  with  J  =  J  +  1,  and  N  others  w.  which  have  that  form  in  I  .   We  wish  to  deter- 
o  1 

mine  the  asymptotic  form  of,  say,  the  v.  in  I  .   To  do  this  we  first  choose  w. 
so  that  its  asymptotic  form  coincides  with  that  of  v.  throughout  I,  up  to  and 
including  terms  of  order  k   .   That  this  is  possible  follows  from  Theorem  I  since 
we  may  choose  the  coefficients  w.  =  v.  for  j  =  0,  ...,J  and  then  determine  w. 
in  I  by  the  construction  of  the  previous  section.   Since  the  component  w. .   is 
obtained  by  integrating  a  first  order  differential  equation,  its  value  can  be 
arbitrarily  prescribed  at  one  point.   Therefore  for  later  convenience  we  prescribe 

(16)     v^(xo)  .  vjf(x0). 

Then  the  theorem  guarantees  the  existence  of  a  solution  w.  with  these  coefficients. 


1 
By  the  linear  independence  of  the  w.  we  may  write 

N 
(17)       v.(x,k)  =   ^   c.s(k)ws(x,k). 


s=o 
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Since  the  coefficients  in  the  asymptotic  form  of  w.  coincide  with  those  of  v. 

through  terms  in  k"  it  follows  that  the  constant  c.   differs  from  the 

Kronecker  5   by  terms  of  order  k    .   Therefore  there  are  constants  7.   such 
is  J  is 

that 

(18)  c.  (k)  =  S.  +7.  ^~J~1+   oCk"*1"1). 

N  1SV        is    IS  N 

When  (18)  and  the  asymptotic  form  of  w  for  x  in  I  are  inserted  into  (17)> 
the  following  asymptotic  form  of  v.  is  obtained  in  I  : 

x  J+l  N         x 

(19)  v.(x,k)  =  exp  k  /  X.(x)dx    )   wS(x)k~S  +  )      exp  k  /   X  (x)dx 

x  s=o  t=l       x 

o  o 


r±t»J(x)  [k"J_1+  oCk-0"-1)]  , 


X  7A4yj_{x)\k.         +  o(.k    )\    ,  x  g  x  . 


Upon  comparing  (19)  with  (15)  we  see  that  whereas  the  asymptotic  form  of  v. 

contains  only  one  of  the  exponentials  in  I  ,  it  contains  all  of  them  in  I  . 

To  determine  c.   we  set  x  =  x  in  (17)  and  solve  the  resulting  linear  equa- 
ls o     x  01 

tions.   They  have  a  unique  solution  because  the  w  are  linearly  independent. 

In  this  way  c.   is  given  in  terms  of  the  v.(x  ,k)  and  w. (x  ,k).   Since  x  be- 
J     is    &  ir   o        1  o  o 

longs  to  both  I  and  I  ,  both  v.  and  w.  have  asymptotic  forms  of  the  type  (15) 
at  x  ,  with  J  =  J  +  1.   By  inserting  them  into  the  expression  for  c  the  asymp- 
totic form  of  c.   is  determined.   In  particular,  to  determine  7..  we  set 

is  '  ' ±t 

x  =  x  in  (19)  and  use  the  asymptotic  form  of  v.  given  by  (15),  which  applies 

at  x  =  x  .   Then  the  coefficients  of  k     in  (19)  yield 

N 

(20)  vJ+1(x   )    =  wJ+1(x    )    +      V       7.+w°(x    )    ,        i    =  1,...,N. 

v  1      v   o  1  o  /_,         it  tx    o      '  ' 

t=l 
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To  solve  (20)  we  recall  that  the  only  nonzero  component  of  w  is  w   because 

o 
v,  is  an  eigenvector  of  the  diagonal  matrix  B  .   Therefore  only  one  7   occurs 
w  o  it 

in  each  component  of  the  right  member  of  (20).   Solving  (20)  yields 

(21)  y.+   =  -ill. 2 ilJL2_  . 

'it  O    ,      X 

vtt(xo} 
From  (l6)  and  (21)  we  see  that 

(22)  7..  =0. 

11 

We  shall  now  compute  the  right  member  of  (21)  in  terms  of  the  coefficients 

B.,  and  thus  express  the  7   in  terms  of  them.   For  this  purpose  it  is  conven- 
J  1  u 

ient  to  introduce  the  special  jump  notation 

(23)  \y{]    =  wJ(xQ+)  -  v^(xQ-)  . 

When  applied  to  B.  or  M.,  the  jump  notation  will  have  its  usual  meaning 

J  J 

(24)  [B.]  =  B.(x  +)  -  B.(x  -)  . 

L  jJ   y   o     y  o 

Let  us  now  consider  (ll)  with  s  =  J  +  1.  We  write  (ll)  for  v.    at  x  - 

x  1  o 

and  for  w         at  x  +  and  subtract 'the  two  equations.      By  hypothesis  M.    and 

BT  .,     .  with  j   =  1.....J  are   continuous.      By  constructidfi  v.(x   )    =  w.  (x   )    for 
J+l-j  °  '        '  J  ivoio 

j   =  0,...,3.      Thus  we  obtain 
J+1t 


(25)  M   [J*]    =    [dvj/dx]    -    [BJ+Jv°   . 


J 

The  left  side  of  (25)  contains  the  jumps  we  need  in  (21)  to  compute  7...  To 
determine  them  we  must  eliminate  the  derivative  from  (25).   To  this  end  we  now 
write  (ll)  with  s  =  J  and  differentiate  it  with  respect  to  x  obtaining 


dM.   T      dv. 

(26)      — i  vJ  +  M.  — i- 
dv   1    1  dx 


,2  J-l    J-l 
d  v. 

l 


dx 


I 


LP, 


dvJ 

l 


-^  vJ  +  BT 

dx    l    J-j  dx 


By  construction  v.(x  )  =  w.(x  )  for  j  =  0,  ...,J  and  by  hypothesis 

dM./dx  is  continuous  if  J  >  0.   Also  by  hypothesis  dB   ./dx  is  continuous  for 
l  J-  j 

j  =  1, ...,J-1  and  B   .  for  j  =  0,...,J-1.   By  the  lemma  dv./dx  is  continuous 

J-J  i 

for  j  =  0,...,J-1.   Thus  upon  subtracting  (26)  at  x  -  from  the  corresponding 


equation  for  w  at  x  +,  we  obtain 

o 


(27) 


M. 
l 


rdvJ1 

i 

_dx  _ 

= 

r,2  J-ll 
d  v. 

l 

- 

_dx  _ 

L  dx2^ 

We  now  multiply  (25)  by  M.  and  combine  it  with  (27)  to  yield 


(28) 


2r  J+ll 


i   L  1  J 


r,2  J-l- 

d  v. 

1 

-( 

[dBjl 

_dx  _ 

L  dx2: 

+  MA+J 


By  continuing  in  the   same  way  we   obtain   for  any  m  §  J  +  1 


(29) 


mr  J+ll 


i*±)m[\]    = 


,m  J-m+1 
d  v. 

1 


dx' 


m 


m-1 


-I 


(M±)" 


r=o 


^-h 


J-m+r+2 


dx' 


m-r-1 


Specializing  (29)  to  m  =  J  +  1  yields 


(30) 


j+i r  J+1t 


(M1)J+1[vi+XJ  = 


r^J+1  o 

a   v. 

1 


cbc 


J+l 


J 

I 


r=o 


(M.)' 


dJ"rB 


r+1 


dx 


J-r 


To  eliminate  the  derivative  of  v.  we  multiply  (30)  by  M.  .   Since  the  only  non- 
zero element  of  v.  is  v ,  . ,  and  since  the  ii  element  of  M.  is  zero,  the  product 
1     ii  1 

,J+1  o/n  J+l 


M  [d   v?/dx   J  =  0.   Then  (50)  becomes 
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J 


(31) 


(M.) 


J+2  r  J+l 


[<+±l    =  '} 


(M.) 


r+1 


rdJ"rB 


r+1 


r=o 


L    dx 


J-r 


x  =  x 


Since  M.    is   diagonal  with  a   zero   in  the   ii  position      (3l)    can  he   solved 
J+l 


for  the   components    [y        |   with  t  /  i.   The  result  is 


(32) 


J+l 


btn--u\-\) 


,r-J-l 


rdJ-rB 


r=o 


r+1, ti 


dx 


J-r 


'ii    ' 


t  t 


Finally,  (32)  and  (21)  yield 


(33) 


7it 


o   J 

T .  . 
11 

O 

tt  r=o 


<V  V 


r-J-1 


rdJ"rB 


r+l,ti 


L   dx 


J-r 


x  =  x   ,  i  /  t. 
o 


This  result  together  with  (22)  completely  determine  7..-   Thus  the  asymptotic 
form  (19)  of  v.  in  I  is  determined.   Then  the  asymptotic  form  of  u . }    the  solu- 
tion of  (l),  is  given  by  multiplying  the  asymptotic  form  of  v.  "by  the  matrix  R, 
according  to  (5) • 

We  may  summarize  our  results  as 

Theorem  II  Suppose  that  in  an  interval  I  of  the  x  axis, 
A(x,k)  has  the  asymptotic  form 


(3*) 


J+l 


A(x,k)  =  V  A.(x)k_,j+  o(k~J_1),     J  g  0 


J=o 


Let  A  (x)  have  N  distinct  eigenvalues  X.(x),  i=l, ...,N, 
o  1 

in  I.   Furthermore,  let  A.  have  J-j  continuous  deriva- 

J 

tives  in  I  for  j=0, ...,J.   Suppose  that  the  only  discon- 
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tinuities  in  I  of  the  derivatives  of  order  (J  +  1  -  j), 
j  =  0, ...,J+1  are  finite  jumps  at  x  .   Let  ui(x)=R(x)vi(x) 
and  uS(x)=RvS(x),  i  =  1,...,N  and  s  =  0, ...,J+1,  where 

vS  and  wS  are  the  vectors  constructed  as  described  above 

i      i 

and  R  is  a  matrix  which  diagonalizes  A  and  has  as  many- 
continuous  derivatives  as  does  A  .  Let  the  constants  y 
be  given  by  (33)  for  i  £   t,  and  7..  =0.   The  coefficients 

uS  and  uS  have  J-s  continuous  derivatives  and  at  x  these 

i      i  ° 

derivatives  of  u.  equal  those  of  u.  for  s  =  0,...,J.   Then 

1  1 

there  exist  N  linearly  independent  solutions  of  (l), 
u  (x,k),i  =  1, ...,N  with  the  asymptotic  forms 


(35) 


u^Xjk) 


exp 


x       -|  rJ+1 
k   /  X  (x)dx     \ 


*-   x 


u  (x)k~  +  o(k    ) 


s=o 


X  S  X 


(36) 


ui(x,k) 


exp 


k   /  Xi(x)dx     ) 


>-   x 


J+l 


u.  (x)k  +  o(k_lJ  J") 


1-  s=o 


)   exp  k   / 


t=l 


X  (x)dx 


^   x 


-0,    >.  -J-l,   /.  -J-ls 
7..u  (x)k    +  o(k    ) 

1  u     0 


X  §  X 


k .        Application 

Let  us  now  consider  the  asymptotic  form  of  the  solutions  of  (2)  when  there 
is  a  finite  jump  in  the  (j+l)-st  derivative  of  the  coefficient  n  (x).   This 
coefficient  is  the  square  of  the  refractive  index  in  wave  propagation.  We  first 
introduce  the  two  component  vector  u  defined  by 


(37) 
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Then  it  follows  from  (2)  that  u  satisfies  (l)  with  A  given  by 

(38)      A(x) 

-n  (x)   0 

We  see  that  A  is  of  the  form  (3)  with  A  =  A,  A .=   0  for  j  11.   The  eigenvalues 

o      j 

of  A  (x)  are  X  =  in(x)  and  X  =  -in(x),  which  are  distinct  if  n  /  0.   If  the 
first  J  derivatives  of  n  are  continuous,  while  the  (j+l)-st  derivative  has  a 
jump  at  x  ,  the  same  is  true  of  A  provided  n(x  )  ^  0.   Thus  the  hypotheses  of 
Theorem  II  are  satisfied  and  therefore  (l)  has  two  linearly  independent  solu- 
tions with  the  asymptotic  forms  (35)  anl  (36). 
If  n(x)  >  0  we  may  interpret  the  factors 


p             X             -1 

X 

ik    /     ndx 

and  exp 

-ik    /      ndx 

u       x            J 

L              X                J 

exp 

o  o 

as  representing  waves  traveling  to  the  right  and  left  respectively.   Then  the 
solution  u_(x)  represents  a  wave  traveling  to  the  left  in  the  region  x  §  x  , 
according  to  (35)-   However,  by  (36)  it  represents  both  left  and  right  travel- 
ing waves  in  the  region  x  g  x  .  We  may  interpret  this  solution  as  corresponding 
to  an  incident  wave  traveling  to  the  left  in  the  region  x  >  x  and  being  partly 
transmitted  and  partly  reflected  at  x  =  x  .   The  reflected  wave  is  of  the  order 
k     with  respect  to  the  incident  wave,  as  we  see  from  (36).   Equation  (36) 
also  shows  that  the  first  component  of  the  reflected  wave,  which  is  kw  according 
"to  (37);  is  proportional  to  7   .   Thus  the  reflection  coefficient  of  the  dis- 
continuity for  incidence  from  the  right  is  proportional  to  7P-,k 

To  determine  7   we  could  use  (33)-   However,  we  shall  instead  analyze  (2) 
directly,  without  transforming  it  into  a  first  order  system  because  the  direct 
analysis  is  usually  simpler  for  a  single  equation.   Therefore,  our  present  cal- 
culations will  be  independent  of  the  preceding  ones.  We  shall  seek  two  linearly 
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independent  solutions  u  (x,k)  and  u  (x,k)  =  u  (x,k),  the  complex  conjugate  of 
u  .   Our  analysis  will  show  that  u_  has  the  asymptotic  form 

x    -t  r  J+l 


(39)      u2(x,k)  =  exp 


-ik   /   ndx 
x 


v.(x)(-ik)"J+  oCk"*1'1) 

(J 


J=o 


x  5  x 


(40)      u2(x,k)  =  exp 


x    i  r  J+l 
/  ndx 
x^    J  Lj=o 


x    -i  r  <j+_l 
f    ndx    V  v.(x)(-ikfJ+  oCk""1"1) 


+  exp 


x 
ik   /  ndx 
x 


7vo(x)(-ik)"J"1+  o(k"J_1) 


X  i  X 


Upon  inserting  (39)  or  (ko)    into  (2)  and  equating  to  zero  the  coefficient  of 
each  power  of  k,  we  obtain 


(*U) 


2n  v!  +  n'v.  =  -v .  , 
J      J     J-l 


The   continuity  of  u     at  x     yields 


j   =  0, ...,J+1 


C*2)  [v  ]    =  0 

J 

C*3>  7  =   -[vJ+J/vo(xo) 


j   «  0,. ..,J 


Continuity  of  u'    at  x     yields 

(kk)  [v]    =   -[v'_J/n(xo)  J«0,...,J 

From  (k2)    and  (Mi-)   we  find 

j   =  0, ...,J-1 
From     (^3)    sad   (^5)   we  obtain 
(^7)  7  -    [vj]/2vo(xo)n(xo) 


(46)  [v1.]    =  0 
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To  calculate  7  we  first  note  that  we  can  find  solutions  of  C+l)  such 
that  v.  has  J  -  j  continuous  derivatives,  for  j  =  0, ...,J.   Then  (^2)  and  (k-6) 
are  satisfied.   Then  we  differentiate  (^l)  m  times  and  subtract  the  result  at 
x  -  from  that  at  x  +  obtaining 

(+8)      2n(xo)[v.W]  +  v.(xo)[n(m)]  =  -  [v.[*+l)],   m=l,  . . . ,  J-  J+l;  j=0,  . . .,  J 

Upon  setting  m  =  J- j+l  in   C+8),   we   obtain  for  j   ^  0, 

(*9)  [v^-J+l)J    =  -    [vj^-j+2)J/2n(xo),      J-a,...,J 

If  we  set  j  =  0  and  m  =  J+l  in  (h-Q)   we  get 


(50)  [vo^]    =   -    [n^]vo(xo)/2n(xo) 


(J+I)j    =  _    Qn(J+l)- 
Now  we   set   j    =  J+l  in    (^9),    and  then  using   (^9)    repeatedly,    we   obtain 


(51)        [vJ+1]  --  [v;]/2n(xo)=     (^ryjc^fj 


/  - 1     \ J+2    /    ^r  (J+Di 


o*    o 


Thus  (^3)  or  (+7)  yields 


On  the  bssis  of  these  results,  we  can  determine  the  reflection  coefficient 
of  the  variable  medium  with  index  n(x).  For  this  purpose  we  must  suppose  that 
n(x)  tends  to  the  value  one  as  x  tends  to  ±  °°.  In  fact,  n(x)  must  tend  to  one 
fast  enough  so  that  n(x)  -  1  is  integrable.  We  now  observe  that  u  represents 
a  wave  incident  from  +  00  which  is  partly  transmitted  and  partly  reflected  at 
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x  =  x  .   The  amplitudes  A.    and  A   „  of  the  incident  and  reflected  waves  are 
o  me      ref 

defined  as 


(53)      A.    =  exp 
w        mc 


k   /  (n-l)dx  -  x 


v0(°°)  > 


w 


r 


A   „  =  exp 

ref 


k   /  (n-l)dx  -  x 


7k    v  (oo) 
o 


Then  the  reflection  coefficient  R  is  given  by 


/   \      ^     rsf    /  .-,  \-J-l 
(55)      Rx  = =  7(-ik)     exp 

inc 


2ik   /  (n-l)dx 


2n  V2ikn(x  ) 


J+l 


-  J+l  - 

d   n 


Ldx 


J+l 


exp 


2ik   /  (n-l)dx  - 


To  determine  the  reflection  coefficient  for  waves  incident  from  the  left, 
we  must  construct  a  solution  which  contains  only  a  transmitted  wave  for  x  >  x 
Such  a  solution  is  u  -  y(-ik)  '   u  ,    at  least  through  terms  of  order  k 
From  it  we  find 


(56)      R2  =  -7(-ik)"  "  exp 


2ik(  /  (n-l)dx  +  x 


2nTx~)  Uikn(x    )]        L   J+l 


J+l  r^  J+l  - 
d       n 


exp 


o 
-2ikM    (n-l)dx  +  xoJ 
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5«    The  case  of  a  discontinuous  A 
o 

Let  us  finally  consider  the  case  in  which  A  is  discontinuous,  vhich  was 

not  included  in  the  preceding  sections.   In  this  case  the  matrix  R  is  not 

necessarily  continuous,  and  therefore  we  do  not  introduce  it.   Instead  we  work 

directly  with  (l)  and  prove  the  following 

Theorem  III  Suppose  that  in  an  interval  I  of  the  x  axis, 
A(x,k)  has  for  large  k  the  asymptotic  form 

(57)      A(x,k)  =  Aq(x)  +  o(l). 

Let  A  (x)  have  N  distinct  eigenvalues  X. (x),  i  =  1,...,N, 

in  I  and  let  A  be  continuous  in  I  except  for  a  finite 
o  * 

jump  at  x  .   Then  there  are  N  linearly  independent  solu- 
tions of  (l),  u  (x,k),  i  =  1, ...,N,  with  the  following 
asymptotic  forms  for  large  k 


(58) 


u.  (x,k)  =  exp  k  /  X.(x)dx  [u.(x)  +  o(l)J, 


x  S  x 

o 


X 

o 


N 


(59)      u.(x,k)  =   ^  exp  k  /\t(x)dx  [c°tu°(x)  +  o(l)],     x  §  xq 

t=l       x 
o 

To  prove  this  result  we  first  apply  theorem  I  with  J  =  0  in  the  region 

x  5  x  to  establish  the  existence  of  N  linearly  independent  solutions  u.  for 

which  (58)  holds.   Then  we  apply  it  again  for  x  g  x  to  obtain  N  other  linearly 

independent  solutions  u.  which  have  there  asymptotic  forms  of  the  type  shown 

in  (58)  with  coefficients  u.(x)  in  place  of  u.(x).   Next  we  express  u.  in  terms 

of  the  u.  in  the  form 

1 
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N 


(60)     u.(x,k)  =   V  cit(k)ut(x,k) 

t=l 


We  then  apply  (60)  at  x  =  x  using  the  asymptotic  form  (58)  for  u  and  the 
corresponding  one  for  u  .   The  resulting  equation  shows  that  c  =  c  +  o(l) 

L-  Xt     It 

where  c   is  determined  by  the  equations. 

IT* 

N 

(61)      u°(xQ)  =   £  cit\(Xo^     ±   "   1*  ••■*■■ 

t=l 

The  solvability  of  (6l)  is  assured  by  the  linear  independence  of  the  u. (x). 

When  the  asymptotic  form  of  c,  .  and  that  of  u, (x)  in  the  region  x  S  x  are 

J   e  it  t  o 

inserted  into  (6o),  (59)  results.  This  completes  the  proof  of  the  theorem, 
vhich  is  essentially  the  extension  of  Theorem  II  to  the  previously  excluded 
case  J  =  -1.  The  results  of  (52),  (55)  and  (56)  of  Section  h  also  hold  for 
J  =  -1  if  2n(x  )  in  them  is  replaced  by  n(x  +)  +  n(x  -). 
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E.  Joppa  Road 

Towson  U,  Maryland 

Attn:  Dr.  D.  M.  Allison,  Jr. 

Director  Engineering  and  Research 

Boeing  Airplane  Company 

Pilotless  Aircraft  Division 

P.O.  Box  3707 

Seattle  2b,  Washington 

Attn:  R.R.  Barber,  Library  Supervisor 

Boeing  Airplane  Company 

Wichita  Division  Engineering  Library 

Wichita  1,  Kansas 

Attn:  Kenneth  C.  Knight,  Librarian 

Boeing  Airplane  Company 

Seattle  Division 

Seattle  1U,  Washington 

Attn:  E.T.  Allen,  Library  Supervisor 

BJorksten  Research  Labs,  Inc. 

P.  0.  Box  265 

Madison,  Wisconsin 

Attn:  Mrs.  Fern  B.  Korsgard 


Convair,   A  Division  of  General  Dynamics 

Corp. 
Fort  Worth,  Texas 

Attn:   K.G.  Brown,  Division  Research 
Librarian 

Convair,  A  Division  of  General  Dynamics 

Corp. 
San  Diego  12,  California 
Attn:  Mrs.  Dora  B.  Burke, 

Engineering  Librarian 

Cornell  Aeronautical  Laboratory,    Inc. 
UllS5  Genesee  Street 
Buffalo  21,   New  York 
Attn:     Librarian 

Dalrao  Victor  Company 
A  Division  of  Textron,    Inc. 
ISIS  Industrial  Way 
Belmont,  California 
Attn:     Mary  Ellen  Addems, 
Technical   Librarian 

Dome  and  Margolin,    Inc. 
29  New  York  Avenue 
Westbury,   Long  Island,   N.  Y. 

Douglas  Aircraft  Company,    Inc. 

P.O.   Box  200 

Long  Beach  1,  California 

Attn:   Engineering  Library  (C-2S0) 

Douglas  Aircraft  Co.,  Inc. 
B27  Lapham  Street 
El  Segundo,  California 
Attn:  Engineering  Library 

Douglas  Aircraft  Company,  Inc. 
3000  Ocean  Park  Boulevard 
Santa  Monica,  California 
Attn:   P.T.  Cline 

Eq.  Sec.  Reference  Files, 

Eq.  Eng.  A250 

Douglas  Aircraft  Company,    Inc. 
2^00  North  Memorial  Drive 
Tulsa,   Oklahoma 
Attn:     Engineering  Library,    D-250 

Electronics  Communication,    Inc. 
1830  York  Road 
Timonium,    Maryland 

Emerson  and  Cuming,    Inc. 
669  Washington  Street 
Canton,   Mass. 
Attn:     Mr.  W.  Cuming 

Emerson  Electric  Mfg.  Co. 

8100  West  Florissant  Avenue 

St.  Louis  21,  Missouri 

Attn:  Mr.  E.R.  Breslin,  Librarian 

Sylvania  Elec.  Prod.  Inc. 
Electronic  Defense  Laboratory 
P.O.  Box  205  -  (Unci) 
Mountain  View,  California 
Attn:   Library 

Falrchild  Aircraft  Division 
Fairchild  Eng.  and  Airplane  Corp. 
Hagerstown,  Maryland 
Attn:   Library 

Farnsworth  Electronics  Company 
3700  East  Pontiac  Street 
Fort  Wayne  lf  Indiana 
Attn:  Technical  Library 

Federal  Telecommunication  Labs. 
500  Washington  Avenue 
Nutley  10,  New  Jersey 
Attn:  Technical  Library 

The  Gabriel  Electronics 
Division  of  the  Gabriel  Company 
135  Crescent  Road 
Needham  Heights  9li,  Mass. 
Attn:  Mr.  Steven  Oalagan 


General  Electric  Advanced  Electronics  Center 
Cornell  University 
Ithaca,  New  York 
Attn:   J.  B.  Travis 

General  Electric  Company 

Electronics  Park 

Syracuse,  New  York 

Attn:  Documents  Library,  B.  Fletcher 

Building  3-l!i3A 

General  Precision  Laboratory,  Inc. 

63  Bedford  Road 

Pleasantville,  New  York 

Attn:   Mrs.  Mary  G.  Herbst,  Librarian 

Goodvear  Aircraft  Corp. 

1210  Masslllon  Road 

Akron  15,  Ohio 

Attn:     Library  D/120  Plant  A 

Granger  Associates 

Electronic  Systems 

966  Commercial  Street 

Palo  Alto,  California 

Attn:   John  V.  N.  Granger,  President 

Grumman  Aircraft  Engineering  Corporation 
Bethpage,  Long  Island,  N.  Y. 
Attn:  Mrs.  A.  M.  Gray,  Librarian 

Engineering  Library,  Plant  No.  $ 

The  Hallicrafters  Company 

UliOl  West  5th  Avenue 

Chicago  2U,  Illinois 

Attn:     LaVeme  LaGioia,    Librarian 

Hoffman  Laboratories,  Inc. 
3761  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Company 
Antenna  Department 
Microwave  Laboratory 
PAilldlng  12,  Room  2617 
Culver  City,  California 
Attn:   M.  D.  Adcock 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 

Attn:  Dr.  L.C.  Van  Atta,  Associate  Director 
Research  Labs. 

Hycon  Eastern,  Inc. 
75  Cambridge  Parkway 
Cambridge,  Mass . 
Attn:  Mrs.  Lois  Seulowitz 
Technical  Librarian 

International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  33,  New  York 

Attn:   Mr.  C.F.  McElwain,  General  Manager 

International  Business  Machines  Corp. 
Military  Products  Division 
Owego,  New  York 

Attn:  Mr.  D.  I.  Marr,  Librarian 
Department  U59 

International  Resistance  Company 
U01  N.  Broad  Street 
Philadelphia  8,  Pa. 
Attn:   Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,  N.  W. 
Washington  7,  D.  C. 
Attn:  Mr.  Delmer  C.  Porta 

Dr.  Henry  Jaslk,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Electromagnetic  Research  Corporation 
711  ldth  Street,  N.  W. 
Washington  5,  D.  C. 
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Lockheed  Aire  reft  Corporation 

2555  N.  Hollywood  Hay 

California  Division  Engineering  Library 

Department  72-75,    Plant  A-l,    Bldg.    63-1 

Burbank,   California 

Attnt      N.  C.   Harnois 

The  Martin  Company 

P.  0.  Box  179 

Denver  1,  Colorado 

Attn:     Mr.   Jack  McCormick 

The  Glenn  L.   Martin  Company 
Baltimore  3,    Maryland 
Attnt      Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,   Maryland 
Attn!     Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,    Rhode  Island 

The  W.  L.  Maxson  Corporation 

U60  West  3lith  Street 

New  York,  N.  Y. 

Attn:     Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,   St.   Louis  3,   Missouri 

Attn:     R.  D.  Detrich,   Engineering  Library 

McMillan  Laboratory,    Inc. 

Brownville  Avenue 

Ipswich,   Massachusetts 

Attn:     Security  Officer,    Document  Room 


Melpar,    Inc. 

3000  Arlington  Boulevard 

Falls  Church,    Virginia 

Attn:     Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,   Massachusetts 
Attn:     N.  Tucker,   General  Manager 
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Pickard  and  Burns,    Inc. 

21:0  Highland  Avenue 

Needhaa  9U,   Mass. 

Attn:     Dr.   J.  T.  DeBettencourt 

Polytechnic  Research  and  Development 

Company,    Inc. 

202  Tillary  Street 

Brooklyn  1,  New  York 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn:  Dr.  John  Ruze 

Radiation,  Inc. 

P.  0.  Drawer  37 

Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.  S.  Carter,  Lab.  Library 

RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attnl  Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldridge  Corporation 
P.O.  Box  1:51:53  Airport  Station 
Los  Angeles  u5,  California 
Attn:  Margaret  C.  Whitnah, 
Chief  Librarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,  Marvland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Library  -  unclassified 

3age  Laboratories 
159  Linden  Street 
Wallesley  81,  Mass. 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Library 

Sandia  Corporation,  Sandia  Base 

P.O.  Box  5800,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Turnbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn:  Library,  Engineering  Division 

Sylvania  Electric  Products,  Inc. 
100  First  Avenue 
Waltham  5b,  Mass. 

Attn:  Charles  A.  Thornhill,  Report  Librarian 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
1U852  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:     M.  L.   Henderson,   Librarian 

A.  S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn:   A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 


Microwave  Radiation  Company  Inc. 

19223  South  Hamilton  Street 

Gardena,  California 

Attn:  Mr.  Morris  J.  Ehrlich,  President 

Chance  Vought  Aircraft,  Ino. 
9311;  West  Jefferson  Street 
Dallas,  Texas 
Attn:  Mr.  H.  S.  White,  Librarian 

Northrop  Aircraft,  Inc. 
Hawthorne,  California 

Attn:  Mr.  E.  A.  Freitas,  Library  DeptJLUS 
1001  E.  Broadway 

Remington  Rand  Univ.  -  Division  of  Sperry 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,  Pennsylvania 
Attn:  Mr.  John  F.  McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,  Inc. 
1221i  Lakewood  Boulevard 
Downey,  California 
Attn:  Engineering  Library  1:95-115 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
Los  Angeles  1:5,  California 
Attn:  Engineering  Technical  File 

Page  Communications  Engineers,  Inc. 
710  Fourteenth  Street,  Northwest 
Washington  5,  D.  C. 
Attn:  Librarian 

Philco  Corporation  Research  Division 

Branch  Library 

1:700  Wissachickon  Avenue 

Philadelphia  Uu,  Pa. 

Attn:  Mrs.  Dorothy  S.  Collins 


Director,  DSAF  Project  RAND 

Via:  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wayland,  Mass. 
Attn:  Mr.  Robert  Borts 

Raytheon  Manufacturing  Company 
Wayland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  0.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  N.  Y. 
Attn:  Engineering  Library 

Thru:  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmingdale,  Long  Island,  N.Y. 

Rheem  Manufacturing  Company 
9236  East  Hall  Soad 
Downey,  California 
Attn:  J.  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  Box  31:6 
Frederick,  Maryland 


Chu  Associates 
p.  0.  Box  387 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Company 
Missile  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50U,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

lbth  Street  and  Constitution  Avenue 

Washington,    D.   C. 

Attn:     Mr.   Otting,    SRY 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  71:6 
Baltimore  3,    Maryland 
Attn:      Engineering  Library 
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Wheeler  Laboratories,  Inc. 
122  Cutter  Mill  Road 
Great  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  nasties  Co. 
Box  91 

Oardena,  California 
Attn:  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
of  the  University  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  li,  California 

Attni  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.  Electronics 
Research  Lab. 

University  of  California 
Electronics  Research  Lab. 
332  Cory  Hall 
Berkeley  U,  California 
Attn:   J.  R.  Whinnery 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
llSOO  Oak  Grove  Drive 
Pasadena,  California 
Attn:  Mr.  I.  E.  Newlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Carnegie  Institute  of  Technology, 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  A.  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Dalman 

University  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:  Prof.  M.  H.  Latour,  Library 

Lihr-ary 

Georgia   Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,    Georgia 

Attn!     Mrs.   J.H.  Crosland,    Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,   303A  Pierce  Hall 

Oxford  Street,   Cambridge  3",   Mass. 

Attn:     Mrs.  E.L.  Hufschmidt,   Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,   Mass. 
Attn:     Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,    Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,    Illinois 

Attn:     Dr.   P.  E.  Moyes,    Department  of 
Electrical  Engineering 

The  Johns  Hopkins  University 
Homewood  Campus 
Department  of   Physics 
Baltimore  18,  Maryland 
Attn:     Dr.  Donald  E.  Kerr 

Sandia  Corporation 
Attn:     Organization  1U23 
Sandia  Base 
Albuquerque,   New  Mexico 


Applied  Physics  Laboratory 
The  Johfls Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:     Mr.   George  L.  Seielstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  20B-221 

Cambridge  39,  Massachusetts 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Document  Room  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J.  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:  Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilanti,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  1U,  Minnesota 

Attn:  Mr.  Robert  H.  Stumra,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E.  Bean 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attnt  Dr.  T.E.  Tice 

Dept.  of  Elec.  Engineering 

The  University  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attnj  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn,  New  York 

Attn:  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn,  New  York 
Attn:  Mr.  A,  E.  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:  Prof.  0.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr.  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:  Prof.  S.  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Morningside  Heights 

New  York,  N.  Y. 

Attn:   Dr.  Schlesinger 

McGill  University 

Montreal,  Canada 

Attn:  Prof.  0.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schulta 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsylvania 

Institute  of  Cooperative  Research 

3b00  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W.  Cumberland  Avenue 
Knoxville  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Att«:  Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  h,    Pennsylvania 
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Ionosphere  Research  Laboratory 

Pennsylvania  State  College 

State  College,  Pennsylvania 

ATN:  Professor  A.  H.  Waynick,  Director 


Institute  of  Mathematical  Sciences 
25  Waverly  Place 
New  York  3,  New  York 
ATTN i  Librarian 
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Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
ATTN:  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.  C. 

ATTN:  Dr.  W.  K.  Saunders 

Applied  Mathematics  and  Statistics  Lab. 

Stanford  University 

Stanford,  (alifornia 

ATTN:  Dr.  Albert  H.  Bowker 

Department  of  Physics  and  Astronomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTN!  Dr.  A.  Leitner 

University  of  Tennessee 
Knoxville,  Tennessee 
ATTN:  Dr.  Fred  A.  Ficken 

Lebanon  Valley  College 
Annville,  Pennsylvania 
ATTN:  Professor  B.H.  Bissinger 

General  Atomic 

P.  0.  Box  608 

San  Diego  12,  California 

Attn:  Dr.  Edward  Oerjuoy 

Department  of  Physics 
Amherst  College 
Amherst,  Mass. 
ATTNi  Dr.  Arnold  Arons 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
ATTN:  Dr.  A.  Erdelyi 

Mathematics  Department 

Stanford  University 
Stanford,  California 
ATTN!  Dr.  Harold  Levine 

University  of  Minnesota 

Minneapolis  l'j,  Minnesota 

ATTNi   Professor  Paul  C.  Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,  California 

ATTN:   Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
ATTN:  Dr.  B.  S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.  C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Mass. 
ATTN:  Dr.  L.  J.  Chu 

Nuclear  Development  Associates,  Inc. 

5  New  Street 

White  Plains,  New  York 

ATTN:  Library 

California  Institute  of  Technology 
Electrical  Engineering 
Pasadena,  California 
ATTNi  Dr.  Zohrab  A.  Kaprielian 


Dr.  Rodman  Doll 
311  W.  Cross  Street 
Ypsilanti,  Michigan 

California  Inst,  of  Technology 
Pasadena,  California 
ATTN:  Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographic  Institute 

Wood3  Hole,  Mass. 

National  Bureau  of  Standards 
Boulder,  Colorado 
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